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Abstract 

Rare b decays provide a unique opportunity to measure Standard Model param- 
eters and probe beyond the Standard Model. We review here the experimental 
progress made in measuring these decays, and the importance of future mea- 
surements, including the possible observation of CP violation. 

I. Introduction 

The dominant decays of the b quark are charged current couplings via a W~ to a c 
quark as shown in Figure 1(a). There are also rare decays to a u quark. Observation 
of these decays has led to measurements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements |V^,| and \V u b\ [1]. 

The b quark can make transitions in other ways. B° — B° mixing, the process 
where a particle changes into its antiparticle, occurs via a "box" diagram with virtual 
W bosons and t quarks inside the box [Figure 1(b)]. The box diagram gives rise to 
large fractions of mixed events: 17% for B° and 50% for B s mesons. 

Flavor- changing neutral currents lead to the transitions 6 — > s and b — s- d. These 
can be described in the Standard Model by one-loop diagrams, known as "penguin" 
diagrams, where a W~ is emitted and reabsorbed [2]. The first such process to be 
observed was b — > 57, described by the diagram in Figure 1(c), where the 7 can be 
radiated from any charged particle line. Another process which is important in rare 
b decays is b — > sg, where g designates a gluon radiated from a quark line [Figure 
1(d)]. A third example of such processes is the transition b — > s£ + £~ which can occur 
through the diagrams shown in Figures 1(e) and 1(f). We consider the loop processes 
shown in Figures l(c)-(f) to be among the most interesting and important rare b 
decays. 

The decay amplitudes for the diagrams shown in Figure 1 are proportional to 
the CKM matrix elements present at each vertex. For the loop diagrams there are 
additional factors of a if a 7 is radiated and a s if a gluon is radiated, as well as a 
kinematic factor which is a function of (m q /mw) 2 - Since the heaviest quark is the top 
quark, it is usually the amplitude involving the top quark that dominates in decays 
via loop diagrams. 
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Figure 1: Feynman Diagrams for b Decays 



B meson decays presently provide the only experimental evidence for penguin 
decays. Although neutral kaons have long been known to mix, penguin contributions 
to kaon decay are hard to identify since the s — >■ d transition leads to the same final 
states as the s — > u transition [3]. Only the dilepton final states can give direct 
evidence for penguins in kaon decay. The channels K + — ¥ ir + vv and K° — ¥ ir°e + e~ 
are the most promising, but the predicted branching ratios are small [4], and thus far 
they have not been observed. 

Loop diagrams in charm decays are suppressed either by CKM matrix elements 
or by small values of (m g /m^) 2 , since the heaviest available quark is the b quark. 
The mixing diagram for D decay is proportional to |K*K/&| 2 X (mb/mw) 2 for b- 
quark exchange and to |K«K/s| 2 X (m s /mff) 2 for s-quark exchange. The same kind 
of suppression factors apply in the case of the penguin diagrams. As a result decays 
such as D — > p7 are expected to be dominated by long distance contributions such 
as rescattering from the related hadronic decay D — > pp. In the case of mixing a long 
distance effect would be D° — > K + K~ — ¥ D°. If these long distance effects are not 
too large, rare charm decays may be sensitive to non-Standard Model effects, since 
the Standard Model predictions for the loop diagrams are so small [5]. 

It is expected that CP violation will be significant in rare b decays. Charge 
conjugation, C, changes a particle to an anti-particle, while parity, P, changes left- 
handed particles to right-handed or vice- versa. In 1964 it was found that the combined 
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operation, CP, showed an asymmetry in neutral kaon decays [6]. CP violation is a 
necessary ingredient in explaining why our local position in the Universe consists of 
matter rather than anti-matter, and thus why we exist. It is of great importance to 
find out whether or not the Standard Model can quantitatively describe CP violation 
in the B system. 

In the parameterization of Wolfenstein [7], the CKM matrix can be described by 
four independent parameters A, A, p and r\. The matrix is given in equation (1). 



Vi 



t l-A 2 /2 
-A 

V AX 3 (l-p-i v ) 



A AX 3 (p - i v ) 
- A 2 /2 AX 2 
-AX 2 1 



(1) 



The A and A parameters have been measured in semileptonic decays of s and b 
quarks. Although r\ and p have not been determined separately, constraints on these 
parameters are given by measurements of the e parameter describing CP violation in 
decay, and by B° — B° mixing and semileptonic b — > Xlv decays. An analysis 
of the allowed parameter space is shown in Figure 2 [8] . Overlaid on the figure is a 

triangle that results from the requirement V^V^. = 0, i.e. that the CKM matrix be 
unitary. Measurements of CP violation in B decays can, in principle, determine each 
of the angles a, (3 and 7 of this triangle independently. 




Figure 2: 
measured 
region is g: 



The CKM triangle overlaid upon constraints in the p — r\ plane, from 
values of V u b/V c b, B° — B° mixing and e in the K° system. The allowed 
iven by the intersection of the three bands. 



In this paper we will review the experimental data on rare b decays and compare 
it with the Standard Model predictions. Following this we will discuss the sensitivity 
of the data to extensions of the Standard Model. Finally we discuss the importance 
of CP violation and the propects for experimental measurements. 
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II. B° - B° Mixing 

The transformation of a B° meson into a B° meson can occur via the diagram 
shown in Figure 1(b). As this is not the main topic of this paper we give only a 
brief summary here and refer the reader to an excellent review for more details [9]. 
The variable that is measured by experiments is x = AM/T, where AM is the mass 
difference between the light and heavy neutral B mesons. The CKM elements are 
related to x via 

X = ^ B Bf 2 B^ B \V t lV td \'F VQCD , (2) 

where Gf is the Fermi constant. The constant Bb and the B meson decay constant 
fs have been calculated theoretically, but the large uncertainties in these calculations 
limit the ability to extract the CKM element \V t d\ from the measurement of a;. To 
determine x experiments either measure the ratio of mixed events to total events inte- 
grated over time (ARGUS and CLEO), or they measure the explicit time dependence 
(ALEPH and OPAL). The extracted x values are shown in Table 1. 

Table 1: x = AM/T Values from B° mixing measurements 



Experiment 


X 


CLEO[10] 


0.65±0.10 


ALEPHfll] 


0.76±0.12 


0PAL[12] 


0.73±0.14 


ARGUS[13] 


0.75±0.15 


AVERAGE 


0.71±0.06 



The band in Figure 2 is derived from equation (2) by assuming Bb = 1, and 
taking an fs range of 160-240 MeV that corresponds to recent theoretical estimates. 
The fraction of mixed events is given by 

X= 2(TT^) (3) 

The measurements of x correspond to a % value of 17% for B° events. There are also 
predictions and experimental limits on mixing in the B s system indicating that the 
mixing has an almost maximal value of 50% [9]. 

III. Observation of Radiative Penguin Decays 

Figure 1(c) depicts the process b — > 57, where the photon can be radiated by 
any charged object in the diagram. This decay is uniquely described in the Standard 
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Model by a "penguin" diagram, with corrections from other diagrams, often called 
"long distance" effects, estimated to be only a few percent (see Section V(A) for 
a detailed discussion). The inclusive process b — > 57 leads to many exclusive final 
states where the s quark hadronizes with the spectator quark. Angular momentum 
conservation forbids the decay B — > K~j, but it is expected that K*(892)^ will be 
a significant fraction of the inclusive rate. The remaining inclusive rate comes from 
higher mass K* resonances and non-resonant K(nir) final states. There are large 
variations among the theoretical predictions for the fraction of b — > 57 that hadronizes 
as B ^ K*f. 

A) Observation of B — >• K*j 

The first successful search for b — > 57 by the CLEO collaboration was for the 
exclusive K**y final state [14]. This is much easier than trying to measure the inclu- 
sive branching ratio for b — > 57, because the final state is completely kinematically 
constrained, and the analysis is similar to that used for reconstructing hadronic B 
meson final states at the T(45') [15]. Neutral clusters in a Csl calorimeter are selected 
with energies between 2.1 and 2.9 GeV, if they have a shower shape consistent with a 
single 7, and if they cannot be combined with another 7 to form a 7r°. The K*(892) 
candidates are searched for in three channels: K*° — > K + ir~, K*~ — > K~ir° and 
K*~ — > K°ir~ . If the energy sum of the K* and the 7 is within 75 MeV of the known 
beam energy Eb eam , then the beam constrained invariant mass 

m B = ^El am - (P^ + K) 2 (4) 

is plotted for each candidate event and an excess is looked for at the known B meson 
mass. 

The difference in shape between jetlike continuum events and spherical BB events 
is exploited by making cuts on several event shape variables to suppress the continuum 
background. The most useful variables are the angle of the thrust axis of the rest 
of the event relative to the candidate thrust axis (cos#r), the second Fox- Wolfram 
moment (-R2) [16], and the sum of the momenta in a 90° cone perpendicular to the 
candidate axis (s_i_) [14, 17]. There is a significant background due to initial state 
radiation (ISR). To suppress this background the events are transformed to the rest 
system of the e + e~ following the emission of the photon. In this primed frame the 
variables cos8' T and R' 2 are recalculated. 

In 1.4/6 -1 of T(45') data there are eight K* ^ and five K*~~f candidates within 
6 MeV of Mb- The continuum background level is one event in each of K*° — > K + ir~ 
and K*~ — > K~ir , and zero in K*~ — > K°ir~ where there are two candidates. This 
is a clear signal for the decay B — > K*^ (Figure 3). The yields of B° — > K* ^ and 
B~ — > K*~~f are consistent. If the relative fractions of B~ and B° produced at the 
T(45) are assumed to be equal, the average branching ratio is (4.5 ± 1.5 ±0.9) X 10~ 5 . 
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Figure 3: Beam constrained mass distribution in GeV for the B — > K*^ candidates, 
dark shaded K*° -> K+ir', light shaded K*~ -> K'tt , unshaded K*~ -> iif 7r- 

B) Measurement of Inclusive b — >• 57 

Recently, the CLEO collaboration has also made the first measurement of the 
inclusive 6 — > 57 branching ratio [18]. The signature for the inclusive process is a 
photon with energy between 2.2 and 2.7 GeV. This region contains 75-90% of the 
signal according to calculations that include the smearing due to the Fermi motion 
of the quarks in the B meson, and the motion of a B meson produced at the T(45'). 

There are large backgrounds to the inclusive signal from continuum jets (e + e~ — > 
qq) and initial state radiation (ISR). These backgrounds are suppressed by two meth- 
ods: a shape variable analysis using a neural network, and a B reconstruction analysis. 
After these cuts have been made, the remaining continuum background is subtracted 
using scaled off-resonance data. There are also small backgrounds from other B de- 
cays, mostly consisting of photons from 7r° and r\ decays that survive a 7r°(7/) mass 
cut because the other photon was not found. As a first approximation these are taken 
from a Monte Carlo simulation, but then a correction is made for any differences that 
are observed between the 7r° and r\ spectra measured in data, and those predicted 
by the Monte Carlo. This takes into account any omissions in the Monte Carlo (e.g. 
b -> sg). 

The neural network analysis uses a set of eight variables defining the event shape. 
The variables i? 2 , s±, R' 2 and cos8' T are as defined in the previous section. In addition 
the energies in 20° and 30° cones parallel and antiparallel to the 7 direction are used. 
The energies in the "away" cones relative to the 7 are found to be particularily useful 
in discriminating against both qq and ISR backgrounds. However since the eight 
variables are highly correlated, and none of them has clear discriminating power 
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compared to the others, they are combined into a joint variable, r, which tends to +1 
for signal, and -1 for continuum background. A neural net is used for this purpose 
since it is the best method of taking into account the correlations between the shape 
variables. 

The B reconstruction analysis combines the high energy photon with a candidate 
X s system, where X s contains either a K s ir + ir~ or a charged track consistent with 
a kaon, and an additional 1-4 pions, of which one may be a 7r°. To be accepted the 
reconstructed decay candidate must satisfy a thrust axis cut, |cos#r| < 0.7, and a % 2 
cut on the combined AE and ms information. If there is more than one candidate 
per event, the one with the smallest % 2 is selected. The reconstruction ambiguities 
usually have the same high energy photon, but different X s systems. This is not 
important if the method is used only to suppress continuum background, and no 
attempt has been made to obtain the corrected X s mass distribution in the CLEO 
analysis. Figure 4 shows the apparent X s mass distribution, with a fit indicating the 
presence of a large component from if* (892). With larger data samples it will be 
possible to study the X s mass distribution and obtain additional information about 
the exclusive decay modes that contribute. 




Figure 4: Apparent Xs mass distribution from the B reconstruction analysis. The 
solid curve is a fit to the expected distribution from a spectator model. The dashed 
curve shows the non-K*(892) component of the fit. 

The two methods for suppressing continuum are complementary. The neural net 
has high efficiency (32%) but modest background suppression, whereas the B re- 
construction method has low efficiency (9%), but suppresses the background by an 
additional factor of 14. According to Monte Carlo studies they should be equally 
sensitive and only slightly correlated with each other. Figure 5 shows the photon en- 
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Figure 5: Photon energy spectra from the neural net analysis, (a) & (b), and from 
the B reconstruction analysis, (c) & (d). In (a) & (c) the on resonance data are 
the solid lines, the scaled off resonance data are the dashed lines, and the sum of 
backgrounds from off resonance data and b — > c Monte Carlo are shown as the square 
points with error bars. In (b) & (d) the backgrounds have been subtracted to show 
the net signal for b — > 57. The solid lines are fits of the signal shape using a spectator 
model prediction. 
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ergy spectra from the two analyses. In Figures 5(b) and (d) the signal shape is taken 
from a spectator model prediction [19]. There is a small excess above the kinematical 
endpoint in Figure 5(b) that is attributed to a statistical fluctuation in the continuum 
background. The measured branching ratios are B(b — > 57) = (1.88 ± 0.74) X 10~ 4 
from the event-shape analysis and B(b — > 57) = (2.75 ± 0.67) X 10~ 4 from the B 
reconstruction analysis. The average result, after taking into account the small cor- 
relations between the two analyses, is B(b — > 57) = (2.32 ± 0.57 ± 0.35) X 10~ 4 , where 
the first error is statistical, and the second systematic. Details of the contributions 
to the systematic error can be found in [20]. 



IV. Theory of Radiative Penguin Decays 
A) Standard Model Prediction for b — >• sj 

The partial decay width for b — > 57 is given by [21, 22]: 

T ( b ^^ = ^^\ V ^ C r ff ^\ 2 (5) 

Since the quark mass, ra&, is not well known, the ml dependence is removed by 
normalizing to the decay rate for b — > civ: 

Y(b — >■ civ) \V cb y 6Trg(m c /m b ) 

where the factor j(m c /m(,) corrects for phase space. In these expressions C^^(fi) is 
an effective coefficient of the electromagnetic loop operator: 

Or = J- mb 3 a a^(l + l5 )b a F^ (7) 

The value of C7 can be calculated perturbatively at the mass scale fi = Mw The 
explicit expression for C~j{Mw) as a function of (m^/M^r) can be found in [23]. The 
evolution from Mw down to a mass scale fi = nib introduces large QCD corrections. 
These are calculated using an operator product expansion based on an effective Hamil- 
tonian: 

n eff (b - , 7 ) = -2V2G F v*v tb £ CiWOiW (8) 

i=l 

Renormalization of the coefficients, (7;, and operator mixing, lead to a value of C^^(fi) 
significantly larger than C~j{Mw) [23]. This increases the predicted rate for b — > 57 
by a factor of 2-3. 

Evidently the prediction for the rate is very sensitive to the QCD corrections. 
The leading log calculation is uncertain to about 25%, primarily because it is unclear 
at which renormalization scale, fi, the effective coefficient, C^^(fi), resulting from 
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the operator product expansion, should be evaluated. Values between fi = |m& 
and fi = 2ra& have been suggested. A next-to-leading order calculation requires the 
evaluation of additional two-loop diagrams, as well as some three-loop diagrams. It 
is hoped that these calculations can be done, since they are expected to reduce the 
uncertainty in the Standard Model prediction to about 10%. 

B) Comparison between b — >• sj Experiment and Theory 

The leading log prediction for B(b -> sj) is (2.8 ± 0.8) x 10" 4 [22, 23]. If 
the next-to-leading order terms that have been calculated are included they tend to 
reduce the prediction to about 1.9 X 10~ 4 [24]. Both these predictions are in excellent 
agreement with the experimental result of (2.3±0.6±0.4) X 10~ 4 . Since the theoretical 
uncertainties are dominated by the choice of the renormalization scale, /i, it is difficult 
to obtain useful constraints on other Standard Model parameters such as m t and V ts - 
The combined CDF and DO measurement of m t = (180 ± 12) GeV [25] is well within 
the range required for consistency with b — > s*y. Ali et al. [26] have set bounds on 
V ts : 

0.62 < ^4 < L1 ° ( 9 ) 

\V c b\ 

This ratio is expected to be one if the CKM matrix is unitary. 



Table 2: Predictions for the ratio of B — > K*^ to b — > 57. 



Author(s) 


Reference 


Method 


B — > if*7 Fraction 


Altomari 


[27] 


Spectator Quark Model 


4.5% 


Deshpande & Trampetic 


[28] 


Relativistic Quark Model 


6 - 14% 


Aliev et al 


[29] 


QCD Sum Rules 


39% 


Ali & Greub 


[19] 


Spectator Quark Model 


(13±3)% 


O'Donnell & Tung 


[30] 


Heavy Quark Symmetry 


10% 


Ball 


[31] 


QCD Sum Rules 


(20±6)% 


Atwood & Soni 


[32] 


Bound State Resonances 


1.6 - 2.5% 


Bernard, Hsieh & Soni 


[33] 


Lattice QCD 


(6.0±1.2±3.4)% 


UKQCD collaboration 


[34] 


Lattice QCD 


15 - 35% 



The fraction of the inclusive b — > 57 rate hadronizing as B — > K*^ depends on 
the B — > K* form factor. This has been calculated by many authors using either 
QCD sum rules, Lattice QCD, or Heavy Quark Effective Theory (HQET). Table 2 
summarizes these predictions for the ratio of B — > K*^ to b — > 57. It can be seen 
that the predictions range from a few percent to 40%. The data suggest a value of 
(21 ± 7)% for this ratio, which is not accurate enough to limit the range of acceptable 
form factor models. It has been suggested by Isgur [35] that the discrepancies between 
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the models could be resolved by using the measured D — > K* form factors as a 
basis for calculating all heavy to light quark form factors. Larger data samples will 
make it possible to distinguish between the predictions in Table 2 and improve our 
understanding of the B — > K* form factor. Until accurate predictions are available 
for the exclusive measurements, they are not as useful as the inclusive measurement 
for constraining the Standard Model or new physics. 

C) Effect of Extensions of the Standard Model on b — >• sj 

The measurement of b — > 57 has inspired a large number of theoretical investigations 
of extensions of the Standard Model that could lead to significant changes in the 
predicted rate for b — > 57. These studies use the upper and lower limits (95% C.L.): 

1.0 x 1(T 4 < B(b^ 57) < 4.2 x 1(T 4 (10) 

to constrain the allowed parameter space of the Standard Model extension being 
considered. Among the most widely discussed models are Higgs doublets, Supersym- 
metry, anomalous WW7 couplings, and anomalous top quark couplings. We give a 
brief summary of these cases below. For investigations into other theoretical ideas 
such as leptoquarks, a fourth generation and left-right symmetric models the reader 
is referred to the review article by Hewett [36]. 

In two-Higgs doublet models there is a charged Higgs that can be inserted into 
the loop instead of the W boson. There are two models for the couplings of the Higgs 
doublets to the quarks, depending on how the fermion masses are generated. In both 
cases the free parameters are the charged Higgs mass, Mjj+, and the ratio of the 
doublet vacuum expectation values, tan/3. With Model I couplings, the b — > 57 rate 
is enhanced at low tan/3, suppressed for values of tan/3 between 0.5 and 1.0, and is 
rather insensitive to large tan/3. Model II couplings always enhance the b — > 57 rate. 
In this case the experimental upper limit requires Mjj+ to be at least 240 GeV even 
for large values of tan/3 [37]. 

Supersymmetry introduces many additional particles that can appear inside the 
loop. In the limit of exact supersymmetry these additional contributions cancel the 
Standard Model contribution, and b — > 57 does not occur at all. In supersymmetric 
models there are charged Higgs bosons with Model II type couplings that enhance the 
rate for b — > 57. Contributions in which down type squarks and either neutralinos or 
gluinos are inserted into the loop are usually found to be negligible. However, there 
are significant contributions when an up type squark and a chargino are inserted 
in the loop. There are several recent analyses of the size and sign of the chargino 
contributions relative to the Standard Model and charged Higgs contributions. It 
appears that there are some regions of the parameter space where the supersymmetric 
model predicts a rate comparable to or below the Standard Model, even for small 
values of M H +. This requires a small stop quark mass, a large value of tan/3, and a 
higgsino mass parameter fi < [38]. 
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D) Constraints on Anomalous Couplings 

The existence of anomalous couplings at the WW^ vertex can be constrained by 
tree-level processes such as e + e~ — > W + W~ and pp — > W~f, and by loop diagrams 
in processes such as b — > 57 [39]. The anomalous couplings are described by two 
parameters, A and Ak, which are zero in the Standard Model, but can acquire non- 
zero values in some extensions of the Standard Model. They modify the value of 
CV(M^), and hence the predicted rate for b — *■ 57: 

C 7 (M W ) = C 7 (M W ) SM + A 1 Ak + A 2 X (11) 

The coefficients Ai and A 2 are functions of (ra^/M^). Since Ai is larger than A 2 for 
m t = 180 GeV, b — > 57 is three times more sensitivity to An than to A. Figure 6 

2.0 

1.0 
< 
-1.0 
-2.0 

-5 -4 -3 -2 -1 1 2 3 4 5 

Ak 

Figure 6: Limits on anomalous WW^ couplings. The shaded regions are allowed by 
the b — > 57 measurement. The region between the shaded regions is excluded by the 
lower limit, the outer unshaded regions by the upper limit. The ellipse shows the 
limit obtained at the Tevatron by the DO experiment (CDF has a similar limit). 

shows the bounds that can be set on A and An from existing data. The limits from 
b — > 57 are complementary to the limits obtained at the Tevatron [40]. Large positive 
and negative values of An and A are excluded by the upper limit on b — > 57. In the 
region around An = — 1 there is a complicated interference between the three terms 
in equation (11). This leads to the exclusion of a narrow band by the lower limit on 
b — > 57. 

Anomalous top quark couplings have also been considered [41]. The first possibil- 
ity is that there are anomalous tt*y couplings in analogy to the WW7 case considered 
above. Once again this would modify C7 through two additional parameters. There is 
also the possibility of anomalous gluon couplings to the top quark that would modify 
Cg, but the constraints on these couplings from b — > 57 are found to be rather weak. 
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Finally there is the interesting point that b — > 57 probes the V-A structure of the 
tbW and tsW couplings [42]. 



V. Searches for Other Radiative Penguin Decays 

A) The Decay B — > pj 

It was suggested by Ali [43] that the ratio of CKM elements could be 

extracted from a measurement of: 




Figure 7: Non-penguin contributions to radiative decays, (a) Color suppressed dia- 
gram for B — > K*(p)^ with a cc intermediate state such as a ip. (b) Color suppressed 
diagram for B° — > p ^ with a uu intermediate state such as a p. (c) Tree level diagram 
for B~ — * p~7 with a uu intermediate state such as a p. (d) Annihilation diagram 
for B~ — > p~7 where the 7 can be radiated from any of the lines. 

Equation (12) is only valid if contributions other than the top-quark loop can be 
neglected in both decay modes. According to Soni [44] there are significant differences 
between the long-distance contributions to b — > 57 and b — > d'y. Examples of such 
additional diagrams are shown in Figure 7. A recent estimate of the long distance 
contributions from virtual ip and p mesons is < 10% for both b — > 57 and 6 — > tZ^ [45]. 
However, there is one contribution from an annihilation diagram (Figure 7(d)), that 
is predicted to be significant by Eilam et al [46], and may be as much as 60% of the 
top-quark loop. Note that this annihilation diagram includes the contribution from 
Figure 7(c) via rearrangement of the quark lines. Since this annihilation diagram only 
applies to B~ decays it is expected that T(B~ — > p~~f) is different from T(B° — > p ^). 
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Deshpande et al. [45] also discuss contributions from the u and c quark loops to 
b — > d'y. These contributions are larger than in b — > 57 and may be as much as 20% of 
the top-quark loop. This would again complicate the extraction of V t d from Equation 
(12). 

CLEO has made a preliminary search for B — > p 7, B° — > p ^ and B° — > u>7 
[47]. A data sample of 2.0/6 -1 at the T(45') results in upper limits between 1.0 and 
2.5 xlO -5 for the three modes. This corresponds to a limit on the ratio in equation 
(12) of 0.34 at 90% confidence level. The search is beginning to be background 
limited. In u>7 the background is primarily from the continuum, whereas in p~~f and 
particularity p°7 there is significant feeddown from misidentified K*^ events. Future 
detectors with better particle identification will be able to suppress this feeddown 
[48], but the continuum background may still be a problem. 

B) Searches for b — >• s£ + £~ 

The process b — > s£ + £~ occurs through a loop diagram with a virtual 7 or Z boson 
(Figure 1(e)), or through a box diagram containing two W bosons (Figure 1(f)). In 
addition the hadronic decays B — > ).fif(*) contribute to the related exclusive decays 
B — > K^£ + £~ through the secondary decays iftO — s- £ + £~. A full understanding of 
b — > s£ + £~ has to include both the short distance contributions from the loop and box 
diagrams, and the long distance contributions from the ip decays, and the interference 
between them [26, 49]. 

At low dilepton masses the dominant contribution from the virtual 7 can be 
directly related to b — > 57. There are also sharp peaks from the ip contributions at 
and m^i which can be directly related to the measurements of the exclusive hadronic 
decays. At high dilepton masses the Z and box contributions are expected to be 
important, as are possible additional contributions from other heavy mass particles. 
The interference between the various diagrams can be studied by measuring the shape 
of the dilepton mass spectrum, and by measuring the lepton-pair asymmetry. 

The high dilepton mass range has been studied at hadron colliders where there is 
a good signature for dimuon pairs. The first search for events with dimuon masses 
between 3.9 and 4.4 GeV was performed by the UA1 experiment [50]. They found 
upper limits of 5.0 X 10~ 5 for the inclusive process b — > sp, + p,~ , and 2.3 X 10~ 5 for the 
exclusive channel B° — > K*°p, + p,~ . Both these limits should be interpreted as referring 
only to the short distance contributions from the loop and box diagrams, since there 
is an extrapolation to the remainder of the phase space under the assumption that the 
long distance contributions are negligible above m^i. Recently the CDF collaboration 
has presented preliminary results from the Tevatron collider. They search over the 
dimuon mass ranges 3.2-3.5 and 3.8-4.4 GeV, and again extrapolate to the full dimuon 
mass range to get upper limits on the short distance contribution of 3.5 X 10~ 5 and 
5.3 xlO -5 for the exclusive channels B° — > K*°p, + p,~ and B~ — > K~ p, + p,~ , respectively 
[51] . 
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In contrast to the hadron collider experiments CLEO has searched for all dilepton 
masses except for the ranges 2.9-3.2 and 3.5-3.8 GeV where the ip( ) contributions 
dominate [52]. The analysis uses standard methods to reconstruct exclusive B meson 
decays from a candidate K or K* meson and a pair of identified leptons. Some typical 
plots of the beam-constrained mass distributions are shown in Figure 8. In an T(45') 
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Figure 8: Beam constrained mass distributions from CLEO for B -> K^l+t: (a) 
K+e+e- (b) K+n+n' (c) K*°e+e 



(d) K* n + n~ 



data sample of 2.0 the background is less than one event in the signal region for 
each of the exclusive channels. The residual background is half from the continuum 
and half from BB events where both B mesons decay semileptonically. Table 2 
summarizes the preliminary upper limits from CLEO for the exclusive channels B — > 
fi + fi and B —> K^e + e . The rate for the decays to electron pairs is predicted 
to be larger than that to muon pairs due to the contribution from low mass pairs 
below the dimuon mass threshold. In some cases the limits from CLEO are close to 
the theoretical expectations. In the future significant increases in statistics at both 
hadron colliders and at T(45") machines are expected to lead to the observation of 
b — > s£ + £~. Eventually there should be enough statistics to measure the dilepton 
mass distribution, and other kinematic variables characterizing the three-body final 
state. Of particular interest is the forward-backward asymmetry of the lepton pair, 
since this is expected to be large in the Standard Model, and may be rather sensitive 
to non-Standard Model physics [26, 54]. 
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Table 3: Results of b — > s£ + £ searches at CLEO. 



B Decay 
Mode 


Candidate 
Events 


Detection 
Efficiency 


90% C.L. 
Upper Limit 


Standard Model 
Prediction[53] 


K+e+e- 


2 


24.4% 


12.0 x 10" 6 


0.6 x 


: 10" 6 


K+n+n' 





15.1% 


9.0 x 10" 6 


0.6 x 


: 10" 6 


K*°e+e- 





9.8% 


16.0 x 10" 6 


5.6 x 


: 10" 6 







5.0% 


31.0 x 10" 6 


2.9 x 


: 10" 6 



VI. Rare Hadronic Decays 

Rare hadronic B decays are described by a combination of a b — > u spectator 
diagram (Figure 1(a)), and a gluonic penguin diagram (Figure 1(d)). Decay modes 
such as 7r + 7r~ and iz^ p T are expected to be described mainly by the spectator diagram 
with a 7r~ or p~ being produced by the W~. A small contribution from a b — > d 
penguin diagram is also expected in these modes. Decay modes such as K~ir + and 
K*~ir + are expected to result mainly from a b — > s penguin diagram, with a small 
contribution coming from the Cabibbo-suppressed spectator diagram where the W~ 
produces a . There are also a few decay modes, such as K°tt° and that are 

described only by a gluonic penguin diagram, and a few modes such as 7r~7r° that are 
described only by a spectator diagram. In modes where both penguin and spectator 
diagrams are significant direct CP violation can occur, as will be discussed in Section 
VII. 

To establish the relative importance of the penguin and spectator amplitudes it 
is necessary to study a large number of decay modes. The first evidence for hadronic 
B meson decays to final states without charmed mesons came from the CLEO obser- 
vation of a signal in the sum of the two decay modes B° — > K + ir~ and B° — > ir + ir~ 
[55]. We will refer to this sum as B° — > h + ir~. Since this publication the CLEO 
data sample has increased by almost a factor of two, and a number of other charmless 
hadronic decay modes have been studied. We also note that there are a few candidate 
events for B° — >■ h + ir~ or B s — >■ h + K~ from the DELPHI and ALEPH experiments 
at LEP [56]. 

A) Decays to and hit Final States 

The signature for a B° decay to two charged tracks is a particularity simple one. 
At the T(4S) the B is almost at rest, and the tracks are back-to-back with momenta 
of about 2.6 GeV. CLEO observed a signal for such events in 1.4 fb~ x of T(4S) data 
[55]. Here we discuss new results from a larger data sample of 2.4 fb~ x [57]. 

The same two kinematical variables are used as in the B — > K*^ analysis (Section 
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111(A)), i.e. the energy sum of the two tracks relative to the beam energy (AE), 
which has an r.m.s. resolution of 25 MeV, and the beam constrained invariant mass 
(ra^), which has an r.m.s. resolution of 2.6 MeV. Particle identification uses dE / dx 
information from the main tracking chamber. Separation between the K + ir~ and 
7r + 7r~ hypotheses comes from the dE / dx information (1.8<r, where a is the rms reso- 
lution), and from the difference in AE (1.7a). The overall separation of 2.5a is rather 
marginal, and is expected to be much better in future detectors [48]. 

The background is due to continuum production of two light quark jets. From 
studies of off-resonance data samples it is known that a cut on the thrust axis, cos#r, 
discussed in Section III(A), is most effective against this background. Requiring 
cos8t < 0.7 removes 95% of the background and only 35% of the signal. There is 
some additional discrimination from the energy distribution of the rest of the event, 
the direction defined by the axis of the two tracks, and the direction of the B meson. 
This information is combined into one variable (J 7 ), using a linear Fisher discriminant 
technique [58]. 

The final signal yields are obtained from a likelihood fit to the four variables, 
AE, tub, T and dE/dx, using an event sample containing the signal region and a 
large sideband in AE and ms from which the background is determined. In the first 
version of this fit all three signal hypotheses are allowed, 7r + 7r~, K + ir~ and K + K~ . 
It is found that the best fit has zero yield for a K + K~ signal. This is expected since 
this decay mode cannot occur via a penguin or spectator diagram. After setting an 
upper limit of 4.0 x 10" 6 (90% C.L.) on B° -> K+K~, a second fit is done in which 
only the first two signal hypotheses are included. The results of this fit projected onto 
the ms and AE axes, are shown as the solid and dotted lines in Figure 9. The event 
histograms result from an event-counting analysis that will be described in section 
VI(B). 

The statistical significance of the signal yield in Figure 9 is determined from the 
probability that the fitted background fluctuates up to the combined yield of signal 
plus background. Although this is determined by Poisson statistics for these small 
event samples, it is conventional to quote the probability in the equivalent number of 
a of a Gaussian distribution. With this definition, the combined significance of the 
two signal modes is quoted as 5<r, with each individual mode having a significance 
of about 2.5<r. These results are interpreted as an observation of the sum of the two 
decays, but not yet as a significant result for either of the individual channels. 

The combined branching ratio for B° — > K + ir~ and B° — > ir + ir~ is measured to 
be (1.8 ± 0.6 ± 0.2) X 10~ 5 . The signal yields and the upper limits on the individ- 
ual branching fractions are given in Table 4. These results are consistent with the 
theoretical predictions given in the last column of Table 4. 

CLEO has made a similar analysis of the decay modes B + — > h + ir°. Here the 
continuum background is larger and the K/ir separation is weaker, since the presence 
of a 7T° leads to a AE resolution of 50 MeV. There are also results for the decay 
modes B° — > 7r°7r°, B + — > K°ir + and B° — > K°ir°. In these three cases only one 
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Figure 9: Projections of the B° — > h + ir~ candidates onto the ms and AE variables. 
The lines show the result of the likelihood fit. In the upper plot the solid line is 
the fitted background, and the dotted line is the fitted signal. In the lower plot the 
lower solid line is the background, the dashed line is the fitted B° — > 7r + 7r~ signal 
the dotted line is the fitted B° — > K + ir~ signal and the upper curve is the sum of all 
three contributions. Shaded events are identitifed as K + ir~, unshaded as ir + ir~. 
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Figure 10: Beam constrained mass distributions for rare hadronic decays to pseu- 
doscalar mesons. The arrows indicate the signal region. 



Table 4: Results of CLEO II Searches for Rare Hadronic B Decays to Two Pseu- 
doscalar Mesons. 

B Decay Signal Yield Fitted B.R. (90% C.L.) Predictions 

Mode (Likelihood Fit) Background xlO" 5 xlO" 5 [22, 59] 



7T + 7r" 



K+ 



7T 



7T + 7T 



K+ir° 



7T 7T 



K°TT + 



K°ir° 



9.4±1? 



7 Q+4-5 
' - y -3.6 



5.ot|:l 



4.9 



+3.6 
-2.8 



1.2tl 



0.9 



C 9+3.5 



o Q + 2.2 
Z -°-1.5 



5.8±0.3 

12.6±0.5 

2.1±0.2 
1.6±0.1 
0.7±0.1 



< 2.0 

< 1.7 

< 1.7 

< 1.4 

< 0.9 

< 4.8 

< 4.0 



1.0-2.6 

1.0- 2.0 
0.6-2.1 
0.3-1.3 

0.03-0.10 

1.1- 1.2 
0.5-0.8 
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signal hypothesis is assumed and no dE / dx information is used in the likelihood fit. 
The mass distributions for these modes are shown in Figure 10, and the results of the 
likelihood fits are summarized in Table 4. It is found that no individual mode has a 
significance greater than 3<r, although most of them are fitted with a small positive 
yield. With the exception of 7r°7r° the other modes in Table 4 are expected to be 
observed with branching fractions comparable to or just below the h + ir~ channels, so 
it is likely that many of these modes will be observed in the near future. 

B) Decays to Vector and Pseudoscalar Mesons 

In this section we discuss searches for the decays B — > irp, B — > Kp, B — > K*ir 
and B — s- K^<f>. The final states for these decays contain three or more particles, 
up to two of which may be 7r°s, and one of which may be a K s . The K* and 
(j) final states are selected by a cut of one natural width about the resonance mass. 
The dE I dx information from the main tracking chamber is used to select the most 
probable decay mode in cases where this is ambiguous. 

In contrast to the previous section, CLEO has used a simple event-counting anal- 
ysis to search for these modes rather than a full likelihood fit. Cuts are made on 
cos#r, F and dE/dx. In addition for decays to a vector and a pseudoscalar meson 
the decay helicity angle, is defined as the angle in the vector meson rest frame 
between the direction of the B meson and one of the decay products of the vector 
meson. Since signal events have a cos 2 ## distribution, a |cos## | > 0.5 cut can be used 
to suppress continuum background. After all these cuts have been made, the event 
yield in the signal region in the ms — Ai? plane is compared to the yield expected 
from an extrapolation of a large two-dimensional sideband region. 

Table 5 summarizes the results from 2.4 fb~ x of T(4S) data. In most cases there 
are few events in the signal region. In the K°p and K^<f> channels there are also few 
events in the sideband region, and we do not quote an estimated background number 
because of the difficulty in extrapolating the yield from such small statistics. There are 
no significant signals in any of the decay modes in Table 5, although 7r ± p :F and K* + tv~ 
do have more events in the signal region than expected purely from background. The 
upper limits on B° — > 7r ± p :F and B + — > K + <f) are close to the theoretical predictions. 

VII. CP violation in Rare Decays 

In the Standard Model CP violation arises from a complex phase in the CKM ma- 
trix, which relates the mass eigenstates to the weak eigenstates. This is an inevitable 
consequence of having three families of quarks. In general, if we have two interfering 
amplitudes, we can write each of them as a product of a strong decay amplitude and 
a weak decay amplitude 
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Table 5: Results of CLEO II Searches for Rare Hadronic B Decays to Final States 
with Vector Mesons. 



B Decay 
Mode 


Event Yield 
(bignal Region) 


Estimated 
Background 


B.R.(90% C.L.) 

^ / 1 (\ — 5 

x 1U 


Prediction 

,,in-5 Too cnl 

xlO [22, 59J 


7T ± / : F 


7 


2.9±0.7 


< 8.8 


1.9-8.8 


7rV° 


1 


1.8±0.6 


< 2.4 


0.07-0.23 


ttV 


4 


2.3±0.3 


< 4.3 


0.0-1.4 


+ 

7T /r 


8 


5.5±1.2 


< 7.7 


iron 

1.5-3.9 


K*+ir- 


3 


0.7±0.2 


< 7.2 


0.1-1.9 


K*°ir° 





1.1±0.3 


< 2.8 


0.3-0.5 


K*+ir° 


4 


1.9±0.7 


< 9.9 


0.1-0.9 


K*°ir+ 


2 


1.0±0.6 


< 4.1 


0.6-0.9 


K+p- 


2 


2.0±0.4 


< 3.5 


0.00-0.20 


K°p° 







< 3.9 


0.01-0.04 


K+p° 


1 


3.8±0.2 


< 1.9 


0.01-0.06 


K°p+ 







< 4.8 


0.00-0.03 


K°<f> 


1 




< 8.8 


0.1-1.3 


K*°<f> 


2 




< 4.3 


0.0-3.1 


K+<f> 







< 1.2 


0.1-1.5 


K* + (j> 


1 




< 7.0 


0.0-3.1 



B = b s e iS >b w e iS ™ . (13) 
Applying the CP operators to these amplitudes results in 

B = b s e i6 °b w e- i6 ™. (14) 

Note that the weak phase has changed sign, while the strong phase has not. The rate 
difference, which may exhibit CP violation is 

r — r = \A + B\ 2 — \A + B\ 2 = 2a s a w b s b w sin(S s — 8 s )sin(S w — 9 W ). (15) 

If two distinct weak decays processes are possible which go via CKM elements with a 
phase difference,then sin(^ — 8 W ) ^ 0. Guaranteeing a strong phase shift, however, 
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is not possible. In fact, the theory of strong decays lacks sufficient power to be able 
to accurately predict the magnitude of such phase differences, or their sign relative 
to the weak phase. 



A) CP violation in B ± Decays to Two Pseudoscalars 

Direct CP violation can occur in charged B meson decays due to interference 
between any two diagrams with different weak and strong phases. In decays to two 
pseudoscalar mesons the penguin and tree diagrams can give rise to integral rate 
asymmetries such as: 

A{K * } " t(b- - + r(5+ - k+*°) (16) 

that are manifestly CP violating. 

There have been several suggestions for measurements of decay rates of B ± mesons 
that could be used to determine the CKM phase sin7 [60, 61]. Although the discussion 
of the derivation of sin7 is complicated, we would like to present the arguments for 
comparing the rates of B ± decays to two pseudoscalar mesons in order to extract sin7 
[60], since it is likely that future experiments will measure these decay rates [48]. 

The final state K~ir° has an amplitude (T s ) from the tree diagram in Figure 
11(a) if the W~ materializes as a K~ , an amplitude (C s ) from the "color suppressed" 
diagram 11(b), and an amplitude (V s ) from the penguin diagram in Figure 11(c). For 
the final state 7r~7r° there are analogous tree (T) and color suppressed (C) amplitudes 
where the s quark in these diagrams is replaced by a d quark. However there is no 
analogous penguin amplitude because the gluon can form dd as well as uu and these 
amplitudes cancel. This is the same as the statement that the I=| 7r + 7r° final state 
cannot be made with a AI=| penguin amplitude. For the K°tv~ final state there is 
only one contribution from the penguin diagram shown in Figure 11(d). 

Assuming SU(3) symmetry, the strange amplitudes are related to the non-strange 
amplitudes by: 

% C s V us fx 



1~ C V u d /,r 

The amplitudes can be summarized as 



0.28 (17) 



A(B~ - tt-tt ) = ~^(T + C) (18) 

A(B~ -> K\-) = V s (19) 

A{B-^K-ir°) = --^{T s +C s + V s ) (20) 

These amplitudes can be related by 

V2A(B~ -> K'tt ) + A(B~ -> K\~) = rV2A(B~ -> t^tt ). (21) 
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Figure 11: Amplitudes contributing to B~ decays to two pseudoscalars. (a) tree 
diagram for K~ir° (b) color-suppressed diagram for K~ir° (c) penguin diagram for 
K~ir° (d) penguin diagram for K°ir~ . The 7r~7r° final state goes through (a) and (b), 
with the W~ — > ud rather than us. 



Amplitude triangles can be constructed from this relationship and its complex con- 
jugate for B + decays (Figure 12). 

Since A(B~ — > K°ir~) involves only one penguin diagram, it is equal to A(B + — > 
K°tt+). On the other hand, 

rV2A(B~ tt-tt ) = are^e-^ (22) 

rV2A(B + tt+tt ) = a r e^ T e i7 , (23) 

where 7 =&Tg(V* b V us ) . Note, that the rates of these two processes are equal since they 
involve a single weak phase and a single strong phase; however, there is a difference 
in phase of 27 between them. In the case of K~ir° and K + ir° the penguin and 
tree contributions interfere and there is a net weak phase shift of arg( V u bV* s V t l V ts )- 

The strong phase shift, 6 = 6q 6p. is also important in determining the actual 

rate asymmetry, A(i^7r°), which has been estimated to be a few percent by several 
authors [62]. CP violation would be explicit if a measured rate difference between 
K~ir° and K + ir° existed, but this requires a strong phase shift as well as the weak 
phase shift. It has been argued that by constructing amplitude triangles as shown in 
Figure 12, the angle 7 can be determined with a twofold ambiguity. Note that if 6 is 
zero 7 can be derived unambiguously even though there is no explicit CP violation 
in B ± -► i^Tr [60]. 

This procedure is only valid if there are no additional contributions other than 
those discussed above. For example the c and u quark penguin loops could gener- 
ate significant asymmetries with different phases from the t quark penguin. There 
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A(KV) 



= A(KV) 



Figure 12: Amplitude triangles relating B ± -> if ± 7r°, B ± -> 7r ± 7r° and B ± -> K ^. 

are several experimental checks that other diagrams can be neglected. The decays 
B~ — > K°ir~ and B + — > K°ir + should have equal rates, as should B~ — > 7r~7r° and 
B + — > 7T + 7r°. If non-equal rates were observed in these decays this would also be an 
observation of direct CP violation, but it would not be possible to extract the angle 
7 from the triangle relation any more. Another check comes from B° — > K + K~ , 
since this cannot be produced by the tree and penguin diagrams of Figure 11. It is 
important that this final state not be observed at a small fraction of the other decays. 

Deshpande and He [63] have pointed out that there is another class of penguin 
diagrams where the gluon in Figure 11(c) and (d) is replaced by a 7, Z, or a box 
diagram. Gluonic penguins couple equally to uu and dd, whereas the 7 and Z pen- 
guins couple differently to u and d quarks. These electroweak penguin contributions 
are expected to be small in B ± — > ^tv , but they could be as large as the tree level 
contributions in B ± — > K ± tv°. If this is the case, the method we have described above 
cannot be used to extract 7 without information from other decay modes. 

There have been two suggestions on how to extract sin7 allowing for possible 
electroweak penguin contributions. Gronau et al. suggest forming an amplitude 
quadrangle including the additional decay mode B s — > rjir [64]. However, measuring 
this rare B s decay appears to be extremely difficult. Deshpande and He suggest in 
a recent preprint the construction of additional amplitude triangles using the decays 
B~ — > K~j]('^ and B + — > K + rj^'^ [65]. The octet part of the rj/rj' system is defined as 

A{K- V8 ) = A{K- V ) cos 9 + A{K- V ') sin 6, (24) 

where 8 is the r\ — r\' mixing angle of about 20° [6] . Two new amplitude triangles can 
be constructed: 

V2A(K-ir°) - 2A(K°ir-) = VQA(K-7] 8 ) (25) 
V2A{K + tt°) - 2A{K\ + ) = V6A{K + V8 ). (26) 
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Figure 13: Amplitude triangles and two solutions a and b for the magnitude of sin7. 
The amplitude A represents A(K°ir~). The amplitude 2 A forms a triangle with the 
amplitudes y/2A(K~ ir°) (line 1) and y/6A(K~r)g) (line 2). This triangle determines 
the amplitude B (line 3). The dotted lines show the equivalent construction for the 
B + decays, leading to the amplitude B (line 4). There is a two fold ambiguity in the 
angle between 3 and 4 due to the possibility of flipping the B + triangle with respect 
to the B~ triangle. This leads to the two results for B — B shown as a and b. 



It is then convenient to construct the amplitude combinations 

B = V2A{K-tt°) - A(K°tt~) (27) 

B = V2A(K+ir°) - A(K\+). (28) 

The complete amplitude construction is shown in Figure 13. The solid lines show 
the amplitude triangle for the B~ decays (equation (25)), while the dotted lines show 
the B + amplitudes (equation (26)). From these triangles B (equation (27)), and 
B (equation (28)), are constructed. The difference B — B (equation (29)), has two 



possible solutions a and b, which are related to the amplitude for B 
to sin7 by: 



7T 7T 



B-B 



i2\/2t 



is 



ud 



A(ir-ir°) 



sin 7, 



and 



(29) 



where 6 is a strong phase shift. From the measured difference from the triangle 
construction a or b, and the measured rate for B~ — > 7r~7r° the angle 7 can be 
determined with a twofold ambiguity. 

The calculation of Deshpande and He assumes that the couplings of kaons, pions 
and etas are related by SU(3), and that the decay amplitudes can be factorized. In 
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equation (29) the amplitude *4("7r~7r°) should be multiplied by a factor (fx I f-w)-, and in 
equations (25) and (26) the amplitude A(K~rj 8 ) by a factor (/„./ f v )- Although there 
is some uncertainty in these theoretical assumptions, the success of factorization in 
explaining B decays to exclusive final states with a D* + and a light hadron has been 
encouraging [66]. 

The angle 7 can also be determined using measured rates in charged B decays 
to D°K final states. The method proposed by Gronau and Wyler [61] uses the three 
related decay modes B~ -> D°K~, B~ -> D°K~, B~ -> DcpK~, where D CP 
indicates that the D° decays into a CP eigenstate, and the corresponding modes for 
B + . The decay B~ — > D°K~ is a Cabibbo suppressed version of B~ — > D°ir~ , while 
the decay B~ — > D°K~ is a color suppressed b — >■ u transition where the virtual 
W~ transforms itself into a cs pair. Interference is possible between these two decays 
modes if the D° decays into a CP eigenstate. Examples of such final states include 
K + K~ , K s ir°i and K s rj. To simplify the discussion only states that are in specific 
angular momentum configurations are used so that their CP is defined as +1 or -1; 
these states are usually denoted as D° and D°. We have 

D ° = 71 l D ° + *1 ■ md D ° = ^\ D °-^\- 

The amplitudes for the three B~ decays modes are related by: 

V2A~(B- -> DlK~) = A(B~ -> D K~) + A(B~ -»■ D°K~). 

Denoting the hadronic phase as S, gives 

V2A~(B- -»■ D^K-) = \A\e i ^' +s) + ~Ae fs . 

The decays to D° need not be equal for B + and B~ , and an asymmetry in them is a 
manifest demonstration of CP violation. A triangle construction serves to determine 
sin7. For more details see [61, 67]. 

B) CP violation in B° — >• 7r + 7r~ due to mixing 

The final state ir + ir~ is one of the simplest in B decay, since there are only two 
pseudoscalar particles in the final state, and there is no spin or angular momentum 
to consider. This final state can be reached from either a B° or a B°, and the s-wave 
production of the two spinless particles means that this is a CP eigenstate. The two 
interfering amplitudes necessary for CP violation are provided by the direct B° decay, 
and the indirect decay following B° — B° mixing. When mixing provides the second 
amplitude for a decay to a CP eigenstate, the strong phase shift disappears from the 
equation relating the measured CP asymmetry to the CKM angles. 

To measure CP violation using mixing we need to make use of the correlated 
production of B°B° pairs. This is done by tagging the number of ir + ir~ events 
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produced with the other B decaying as a B°, as opposed to a B°. An example 
of a suitable tag is lepton flavor. These numbers are time dependent functions of 
T = t — t' (in units of mean lifetime), where t is the decay time of the ir + ir~ and t' is 
the decay time of the other B: 

1Z{T) oc e~' T '(l — sin —^—T sin 2a) (for tt+tt", B°), (30) 

TZ(T) oc e~' T '(l + sin ^y^T sin 2a) (for tt+tt", B°). (31) 
It can be shown that 1Z{T) = 1Z(—T). The CP asymmetry is 

, , , , K(T)-K(T) . Am m . n 

At(tt + tt~) = — — =7—7 = sin T sin 2a (32) 

y J n{T) + n{T) r v ; 

We can measure a time independent asymmetry by integrating over t—t'. However, 
if the C parity of the intial B°B° pair is -1, as is the case in T(45) decay, the time 
integrated rate is zero. This is not the case in hadron colliders, or if the initial state 
is B°*B°. 

Measuring CP violation in the ir + ir~ decay determines sin2a, where 

a = Mg(V ud V: b /V td V t l). (33) 

However, there is a problem due to the presence of a decay amplitude related to 
the penguin diagram shown in Figure 1(d). A ir + ir~ final state is produced when 
the t quark in the loop couples to a d quark rather than an s quark. Although this 
is suppressed it could turn out to have a significant effect on the CP asymmetry 
measurement. Gronau and London [68] have shown how this contribution can be 
isolated by measuring the rates for the 7r + 7r° and 7r°7r° processes. Equations (30) 
and (31) are modified as follows: 

K ( T) oc (l±f. + izJf cos + lm( sin ^|T|) (34) 

R(T) K e-'rl(i^-i^ CO s^T-I™, sin^|T|), (35) 

where £ is a parameter related to a. If only the direct and the mixing amplitude 
are present |£| = 1, and we are left with equations (30) and (31). We can look for 
the presence of the cosine term experimentally by forming a new time dependent 
asymmetry: 



A i + -x ^(|r|) + ^(-|r|)-^(|r|)-^(-|r|) i-|£| 2 Am 

Am[ * * >-^\T\)+n-\T\)+n\T\)+n-\T\)~T+W ~ l l ' ( J 
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The asymmetry A\t\ leads to a non- vanishing asymmetry, A, in the time integrated 
rates, even at the T(45'): 



^ + -> = 'ttw^Wt (37) 

If it turns out that the 7r°7r° rate is comparable to 7r + 7r~, the penguin amplitude is 
important but it should be possible to measure it. If it is small enough to be difficult 
to measure, then the penguin amplitude is likely to be unimportant and a comes out 
simply. 

Another method for extracting the size of the penguin-tree interference is to com- 
pare the integral rates for B° — > 7r + 7r~ and B° — > K + ir~ with those for B° — > 7r~7r + 
and B° — > K~ir + [69]. A recent preprint by Deshpande and He [70] shows that the 
assumption of SU(3) symmetry leads to the following simple relationship between the 
time integrated rate asymmetries for these decays: 

A(7r+7r-) « -|f-A(#+7r-) (38) 
Jk 

and that this result can be used to correct the measurements of the time-dependent 
CP asymmetry in ir + ir~ for the effect of the penguin amplitude. 

C) CP violation in radiative penguin decays 

A measurement of the ratio of b — > d'y to b — > 57 determines |Vtd|/|Vt s | if the 
amplitudes are described by the t quark loop. As discussed in Section V(A) there 
may be other amplitudes that are significant in b — > d'y. If this is the case then these 
additional amplitudes can give rise to direct CP violation in these decays. 

One source of CP violation is the presence of three different loop diagrams involv- 
ing the ii, c or t quarks. Gluon exchange provides the necessary strong phase shifts 
between these loop diagrams. Naively, one would expect that the u and c diagrams 
would be highly suppressed due to the relatively small quark masses. However both 
Deshpande et al. [45] and Soares [71] find significant u and c loop contributions 
in b — > d'y. Soares has explictly calculated the amount of CP violation from these 
sources and finds asymmetries of (2 — 8) X 10~ 3 for b — > 57 and a significantly higher 
value of (2 — 30) X 10~ 2 for b — > d'y. For the exclusive decay modes Greub et al. [72] 
estimate that the CP asymmetries are 1% for B — > if*(890)7 and 15% for B — > p~f. 

There are other diagrams that can provide a source of CP violation. Of particular 
interest is the suggestion that non Standard Model contributions to b — > 57 can lead 
to large CP asymmetries, since the asymmetry predicted by the Standard Model is 
rather small. As an example of such possibilities we mention the paper by Wolfenstein 
and Wu [73] that calculates the expected level of CP violation in a two Higgs doublet 
model. Here the charged Higgs is present in the loop instead of the W~ , and there 
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is an arbitrary phase factor associated with the two Higgs doublets. Depending on 
this phase the asymmetry in b — > 57 could be anywhere in the range 0-10%. If a CP 
asymmetry larger than 1% were observed in b — >■ 57, it would provide strong evidence 
for physics beyond the Standard Model. 

VIII. Conclusions 

Loop diagrams were first discovered in the mixing amplitude for neutral kaon 
decays. Mixing has also been observed in the neutral B mesons [Figure 1(b)]. CLEO 
has established the existence of the radiative penguin decay b — > 57 [Figure 1(c)], by 
measuring the exclusive branching fraction for B — > K*j to be (4.5 ± 1.5 ±0.9) X 10~ 5 , 
and the inclusive branching fraction for b — > 57 to be (2.3 ± 0.6 ± 0.4) X 10~ 4 . These 
values are in agreement with the Standard Model prediction, and set bounds on the 
parameters of some extensions of the Standard Model. Other penguin decays such as 
B ~^ PI an d b — > s£ + £~ have been searched for but have not yet been seen. 

There is also a CLEO measurement of (1.8 ± 0.6 ± 0.2) X 10~ 5 for the branching 
fraction of the sum of the decays B° — > K + ir~ and B° — > ir + ir~. While the K + ir~ 
final state is thought to occur primarily through the gluonic penguin diagram [Figure 
1(d)], the 7r + 7r~ final state occurs primarily via ali-»ii tree level transition [Figure 
1(a)]. The data favor equal branching ratios for the two decay modes, but only 
exclude either one being zero at a level of significance equivalent to about 2.5<r. 

Decays of B mesons to two pseudoscalar mesons are particularly important for 
the study of CP violation. The ir + ir~ final state can be used to measure the angle a 
in the CKM triangle. Corrections for a penguin contribution to B° — > ir + ir~ can be 
made by making additional measurements of B + — > 7r + 7r° and B° — > 7r°7r° and using 
isospin, or by measuring the rate asymmetry between B° — > K + ir~ and B° — > K~ir + 
and using SU(3) symmetry. It is likely that the CKM angle 7 will be measured using 
charged B decays to K^tt , tt^tt , K ^ and K^rj^ \ or using charged B decays to 
D°K ± final states. 

We are looking forward to the measurement of additional rare b decays such as 
B ~^ PI an d b — > s£ + £~ . The Standard Model is already constrained by the existing 
measurements of B° mixing and b — > 57. If deviations from the Standard Model are 
found in other radiative penguin decays, in CP asymmetry measurements, or in K or 
D meson decays, then possible extensions of the Standard Model must both explain 
the observed deviations and be consistent with the other measurements. 

There are planned improvements to the CESR/CLEO symmetric 5-factory, and 
asymmetric B factories are under construction at KEK and SLAC [48], all with pro- 
jected luminosities about ten times higher than that currently achieved at CESR/CLEO 
In addition hadron collider B experiments are being pursued (HERA-B, Tevatron, 
LHC)[74]. We hope that these efforts will lead to the observation of CP violation in 
the B system, and that some evidence for non-Standard Model effects will be found. 
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